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Abstract. We show that the partition functions which enumerate Donaldson-Thomas invari- 
ants of local toric Calabi-Yau threefolds without compact divisors can be expressed in terms 
of specializations of the Schur measure. We also discuss the relevance of the Hall-Littlewood 
and Jack measures in the context of BPS state counting and study the partition functions 
at arbitrary points of the Kahler moduli space. This rewriting in terms of symmetric func- 
tions leads to a unitary one-matrix model representation for Donaldson- Thomas theory. We 
describe explicitly how this result is related to the unitary matrix model description of Chern- 
Simons gauge theory. This representation is used to show that the generating functions for 
f^ I Donaldson-Thomas invariants are related to tau-functions of the integrable Toda and Toeplitz 

Cn ■ lattice hierarchies. The matrix model also leads to an interpretation of Donaldson- Thomas 

^ theory in terms of non-intersecting paths in the lock-step model of vicious walkers. We further 

■ show that these generating functions can be interpreted as normalization constants of a corner 

^_. ^ growth/last-passage stochastic model. 



ji. . 1. Introduction and summary of results 

(— I ' Donaldson-Thomas theory computes enumerative invariants associated to the number of 
points in the moduh spaces of ideal sheaves with trivial determinant on a three-dimensional 

^SJ . Calabi-Yau variety [HE]- The partition functions can be rephrased in terms of the counting of 

^ I noncommutative U{1) instantons in a six-dimensional topological gauge theory OUj. In this 

tH- I way the Donaldson-Thomas partition functions may be regarded as generating functions which 

^O ' count BPS bound states of DO and D2 branes in a single D6-brane, at least for appropriate val- 
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ues of the 5-field. They can also be interpreted combinatorially in terms of the enumeration of 
plane partitions (three-dimensional Young tableaux) with boundary conditions along the three 

f^ I axes given by ordinary partitions (Young diagrams), where the plane partitions are glued to- 

gether along common boundaries to form a crystal configuration. This leads to an interpretation 
of Donaldson-Thomas theory in terms of the statistical mechanics of crystal melting [5]. The 

S^ . melting crystal formulation connects Donaldson-Thomas theory with topological string theory 

through the formalism of the topological vertex [6|. For local toric backgrounds, the generating 
functions for Donaldson-Thomas and Gromov-Witten invariants are related by a simple change 
of variables [3 [8] . 

In this paper we will construct some new statistical mechanical models of Donaldson-Thomas 
theory for local toric Calabi-Yau backgrounds which have no compact divisors. We build upon 
the representation of these partition functions in terms of random Young diagrams, which follows 
implicitly from the original expansion of the topological vertex [6j. To handle random Young 
diagrams, one needs to define a probability measure on the set of Young tableaux. The classic 
example is the Plancherel measure, introduced in the 1970's by Kerov and Vershik [9]. It is of 
interest in field theory and string theory, wherein the partition functions of AA = 2 supersym- 
metric gauge theory in four dimensions |10^ \TT\ and of topological string theory on local toric 
curves [S] can be conveniently expressed in terms of Plancherel and g-Plancherel measures. 



Date: May 2010 HWM-10-22 , EMPG-10-10. 
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More generally, one can define the Schur measure A^schuD introduced by Okounkov [12J, such 
that M-Schuri^} for a partition A is proportional to Sx{x) 5\{y). Here 5\{x) are the Schur 
polynomials, and x and y are two independent (possibly infinite) sets of variables. There are 
several important properties satisfied by the Schur measureo For instance, as we shall see in 
Section [21 after proper specification of the variables it contains the Plancherel and g-Plancherel 
measures [131 [T^ as particular cases. Moreover, the Schur measure has correlation functions 
of determinantal type, which is a common feature of various problems in statistical mechanics, 
enumerative combinatorics and probability theory that leads to their explicit solution. One of 
the remarkable mathematical results of the last decade has been the expression of probability 
measures on partitions in terms of determinantal point processes |15t [T6| I17j . which are very 
often of random matrix type |18tll9j. Some straightforward generalizations of the Schur measure 
that we shall discuss in Section [2] are defined in terms of the Hall-Littlewood polynomials, or 
the Jack polynomials which are one-parameter extensions of the Schur polynomials 5x (x) \20\n 

In this paper our interest in the Schur measure and its generalizations rests in the property 
that the partition functions of Donaldson-Thomas theory, and also of topological string theory 
and of the BPS state counting for D6-D2-D0 branes [21], on the backgrounds we consider admit 
very natural expansions in terms of these measures. More precisely, as we describe in detail in 
Section [21 they are normalization constants of particular cases of such measures. These special 
instances come from the specification of the two independent (infinite) sets of parameters x and 
y in terms of the variables of the generating functions for the Donaldson-Thomas invariants [21] . 
Several mathematical properties of these measures are utilized throughout this paper. Specific 
examples of such rewritings in terms of symmetric functions can be found in [22^ [23] in the 
context of the melting crystal model with external potentials, and in \2^ in the context of 
topological string theory (see also [21] for a review). 

This representation of the partition functions in terms of symmetric functions is especially 
interesting due to the well-known combinatorics theorem of Gessel [25], which shows that the 
normalization constant of the Schur measure can be written as a Toeplitz determinant. By 
exploiting the classical Heine-Szego identity which links Toeplitz determinants with unitary 
matrix models |26j , it then leads to a unitary one- matrix model representation for the Donaldson- 
Thomas partition functions. This is the subject of Section [3l These matrix models appear to be 
different from the ones previously found for Donaldson-Thomas theory [27] and for topological 
string theory [2BJ , which involve infinite-dimensional hermitean multi-matrix integrals with non- 
polynomial potentials. Since our unitary matrix model is characterized by an infinite number 
of eigenvalues and its weight function is a Jacobi theta-function, we can immediately relate it 
to the unitary Chern-Simons matrix models [22 [301 [S] with gauge group C/(oo). This is done 
by completing the unitary matrix model representation of Chern-Simons theory by giving the 
partition function in terms of a Toeplitz determinant. 

The equivalent representations in terms of matrix models, Toeplitz determinants and Schur 
measures imply, among other things, that the generating functions for Donaldson-Thomas in- 
variants are integrable. We study this problem in detail in Section [H Using the work of Sato [32j 
and Segal- Wilson [33] , which show that tau- functions of integrable hierarchies admit expansions 
in terms of Schur functions, we relate the Donaldson-Thomas partition functions to tau- functions 
of the integrable Toda lattice hierarchy. The equations of the hierarchy, together with the string 
and divisor equations, uniquely determine the entire theory. We consider various points of view 
on this issue in gauge theory and string theory which suggest that, in appropriate instances. 



See [12] for its more representation theoretic properties. 

The Macdonald polynomials [20] are the most general symmetric polynomials known and they include all the 
other ones as special limiting cases of their two parameters, but we shall not need them here. 
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the pertinent tau-function is either of 2-Toda or 1-Toda type. For example, in the free fermion 
formulation |221 [23] there is a hidden symmetry yielding reduction to the 1-Toda lattice hier- 
archy, which is related to the integrability structure of super symmetric gauge theory in four 
and five dimensions. On the other hand, the 1-Toda structure hidden in the 2-Toda formal- 
ism is also evident in both ordinary and g-deformed two-dimensional Yang-Mills theory. The 
results of [M] and |35l §1.4.2] seem to hint that double Hurwitz numbers may be more natural 
or better suited to describe the branched cover interpretation of these gauge theories, which 
follows from the philosophy of our treatment as well. This extends the usual description of two- 
dimensional Yang-Mills theory in terms of simple Hurwitz numbers. Insofar as the g-deformed 
gauge theory serves as a nonperturbative completion of the A-model topological string theory 
on certain backgrounds, this integrability structure seems to be related to the fact that equivari- 
ant Gromov-Witten theory can be either of 1-Toda or of 2-Toda type. For example, topological 
string theory on the resolved conifold belongs to the 1-Toda lattice hierarchy, while on toric Fano 
surfaces it belongs to the 2-Toda lattice hierarchy |36j . The latter geometries contain compact 
divisors and so fall out of the class of backgrounds that our main line of development applies 
to in this paper. This explanation for the absence of symmetry reduction of the Toda lattice 
hierarchy appears to be a generic feature of topological string partition functionsjj 

A rather generic problem of determinantal type is that of non-intersecting paths or non- 
intersecting Brownian motion in the continuum. The formalism developed by Karlin and Mc- 
Gregor in the 1950's [39j provided elegant determinantal expressions to describe non-intersecting 
Brownian motion. These formulas were used in [401 HT] to express Wilson loop observables in 
Chern-Simons gauge theory as probabilities in a model of N non-intersecting Brownian motion 
particles, where A'^ is the rank of the gauge group. The system of non-intersecting Brownian mo- 
tion paths was introduced by de Gennes to study chains of polymers under steric constraints |42j . 
Later on, Fisher |43l \S^ introduced two models of non-intersecting (vicious) random walkers 
in order to model domain walls in two-dimensional lattice systems: the lock-step model and 
the random-turns model. The latter model is intimately related to unitary random matrices 
and the Gross- Witten model [40j . It was also shown in [40 [ [IT] how some of the vicious walker 
expressions given by Fisher are related to Chern-Simons observables. For example, the proba- 
bility of reunion for A^ vicious walkers on a line gives the partition function of Chern-Simons 
gauge theory on S^ with gauge group U{N). The problem can be equivalently understood as 
the Brownian motion of a single particle on the Weyl chamber of the gauge group. In |40] . the 
eventual role of the lock-step model in gauge theory was left as an open question. 

The lock-step model of vicious walkers on a one-dimensional lattice allows each walker at the 
tick of a clock to move either one lattice site to the left or one lattice site to the right, with 
the restriction that no two walkers may arrive at the same lattice site or pass one another. 
In Section we will show that our matrix model expressions can be interpreted in terms of 
non-intersecting paths in this model with infinitely many vicious walkers. This follows directly 
from [15] which gives a vicious walker interpretation of the generic matrix averages that describe 
the Donaldson-Thomas partition functions. 



This description appears to be equivalent to the topological sigma-model description of Gromov-Witten 
theory on P^ [37], which is equivalent to a large A*' hermitean one-matrix model with non-polynomial potential 
that belongs to the 1-Toda lattice hierarchy. In contrast, the equivariant Gromov-Witten theory of P^ is governed 
by the 2-Toda lattice hierarchy [38] . 
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In addition, it turns out that the Schur measure employed here is the basis of a generahzed 
version of the corner growth model [18t I19j . and hence we can also interpret the Donaldson- 
Thomas partition functions as normalization constants of this stochastic processo This model 
is intimately related to other random models, such as the discrete polynuclear growth model, 
non-intersecting paths, and random tilings of Aztec diamonds |18t [T9j (see |15j for a review). 
It is believed to belong to the Kardar-Parisi-Zhang universality class for stochastic growth 
processes [15] . 

These two descriptions imply that Donaldson-Thomas theory has an interpretation in terms of 
both discrete and continuous random systems. This gives an indication that the non-intersecting 
paths system of ^45] is intimately related to the continuous last passage model in [19j. Our treat- 
ment of the BPS bound state partition function holds at all points in the Kahler moduli space, 
and hence naturally incorporates its discontinuity due to wall-crossing phenomena. In these 
statistical mechanics models, the jumps across walls of marginal stability have a natural inter- 
pretation in terms of the creation or destruction of particles and independent random variables. 
Furthermore, using results of [13] we shall also see that the Donaldson-Thomas partition func- 
tions are generating functions of certain types of infinite integer matrices that satisfy specific 
symmetry conditions. It would be interesting to relate more directly this enumerative interpre- 
tation to the actual integer Donaldson-Thomas invariants that count ideal sheaves (equivalently 
instantons or D6-D2-D0 bound states). 



Note added. While we were completing the present paper, the preprint [49j appeared, whose 
results overlap with ours, but with a different approach and theme. Their derivation of the 
unitary matrix models is based on observing that the Chern-Simons matrix model is related to 
the MacMahon function by expansion of the weight function of the matrix model. Based on the 
melting crystal formalism, both a free fermion formulation and the Gessel-Viennot determinantal 
expression which counts non-intersecting paths is then employed to find explicit expressions for 
the matrix models. As summarised above, our approach is different and yields these same 
results within a different setting. For example, the non-uniqueness of the weight function of the 
matrix model, discussed after eq. (4.20) in [l9], is emphasised below in ()3.18p and the discussion 
afterwards in a completely different way; this has also been noticed previously for Chern-Simons 
matrix models. Notice that the Lindstrom-Gessel-Viennot theorem used in [39] is essentially 
the same as the Karlin-McGregor theorem that was used previously in [401 [3T] to obtain Chern- 
Simons observables. Although the Karlin-McGregor theorem leads to a two-matrix model, it is 
explained in [31] how to relate it to the one-matrix model of Chern-Simons gauge theory. This 
result is also used by Johansson to find the Schur measure (j5.5p in the study of the generalized 
corner growth model that we employ in Section [5l 

Acknowledgments. This work was supported in part by grant ST/G000514/1 "String Theory 
Scotland" from the UK Science and Technology Facilities Council. 

2. Donaldson-Thomas theory and symmetric function measures on partitions 

2.1. AfRne space. The partition function of Donaldson-Thomas theory in the simplest case of 
the Calabi-Yau threefold C^ is given by the MacMahon function 

oo ^ 

(2.1) Z^^{q) = l[- -^^M{q), 



In [461 147) similar partition functions in two and three dimensions are related to growth processes described 
by the integrable XXZ spin chain and to a generalization thereof. 
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where we choose a string couphng constant gs such that the quantum parameter q = e "^^ 
satisfies \q\ < 1. This is the generating function for plane partitions |50t I51j. They correspond 
to pointhke instantons, or DO-branes, in the D6-brane gauge theory on C^, with charge equal to 
the number of boxes in the plane partition. Each plane partition is equivalent to a monomial 
ideal which corresponds to an ideal sheaf on C^. 

As pointed out by [5l[22l[23l[2T], this partition function can be rewritten in terms of Schur func- 
tions. This expression is the crux of the equivalence between Donaldson-Thomas and Gromov- 
Witten theories in the toric case [71[8l[3], with the latter computed from the topological vertex 
formalism for topological string theory [6] within the melting crystal formulation [5J. However, 
it also follows directly from the enumerative expression. The Schur polynomials sx{x) [51, .20] 
constitute a basis of symmetric functions in a given set of variables x = (xj)j>i and are in- 
dexed by Young diagrams (ordinary partitions) A = (Aj)j>i, with Aj > Aj+i > giving the 
length of the i-th row. If the variables x are regarded as eigenvalues of some matrix M S sin, 
then Sa(x) = Trx{M) is the trace of M in the irreducible s^n-representation associated to A. 
The Schur polynomials may also be more explicitly defined in terms of the skew-symmetric 
polynomials o^ = detij{x^^ ) as Sa(x) = aA(a;)/ao(x). By using the Cauchy identity [51] 

(2.2) J]5.(x)s.(y)=n r:^"^' 

and considering the case of an infinite number of variables with the specializations Xj = g*~^' ^ 
and yj = q^~^''^, then ()2.2p directly gives the expression |21] for the Donaldson-Thomas partition 
function of C^ in terms of Schur functions, 

(2.3) zE^iq) = ^.,{q^-y^f. 

A 

Let us phrase this result in terms of the Schur measure introduced by Okounkov in [12] . 
It assigns to each partition A the weight A^schur{-^} = ^5\{x)5\{y), where 5x{x) are Schur 
functions. Then 

(2.4) Vn{X} = ^ Y. SA(x)sA(y) 

A:Ai<Ar 

is the probability that the number of boxes in the first row of the associated Young diagram 
is < A^. Taking the limit A^ — )• oo, the normalization constant Z is given by (|2.2p . Thus the 
Donaldson-Thomas partition function ()2.3p is given by the normalization constant of the Schur 
measure, specialized at Xi = g*~^'^ and yj = q^^^'"^. 

By making the specification 5^(1, 1, • • • ) = dim A, the Schur measure contains the Plancherel 
measure 

/n rr\ A^ rxT ^ dim A \ 2 

(2.5) A^PlanchjA} = 



|A|! 

where the dimension of the corresponding irreducible representation of the symmetric group 
S\x\, with |A| := ^,- Aj the weight of the representation, is given by 

(2.6) *-^ = nsR- 

MSA ^ ' 

Here h{u) = Xi + X', — i — j + 1 is the hook length of the box u = {i,j) in A, and A' denotes 
the conjugate partition to Ao Our specification is instead q~l'^l/2 s;^(l, q, . . . , g"~^), and we know 



The conjugate of a partition is obtained by interchanging rows and columns in its Ferrers graph. 
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(2.7) SA(l,g,...,g"-^) = g"(") n 



that [SD ED] 

^r^-l^, ^ „n(A) TT [^ + C (u)] 

J- J- [/i(«)] ' 

where n(A) = X]i>i (^ ~ 1) -^« ^^d, for each box u = {i,j) of the diagram A, c(u) = j — i is the 
content of u. The square brackets here denote (;-numbers, 

a/2 _ -a/2 

and the right-hand side of ()2.7p is the g-deformation of the dimensions of sin representations, 
i.e. the quantum dimensions dinig A [52j. Then the Schur measure with this speciaHzation is a 
g-Plancherel measure, essentially the g-deformed Plancherel measure discussed in [T3]|j Thus 
the normalization constant of the g-Plancherel measure gives the Donaldson-Thomas partition 
function on C^. 

2.2. Conifold. There are two non-isomorphic crepant resolutions of the conifold singularity 
zi Z2 — z^Z4 = in C^ [53]. The topological string partition function of the resolved conifold X = 
Opi (—1) © C'pi(— 1) can be expressed in terms of the Donaldson-Thomas generating function 

(2.9) Z^.,{q,Q)=M{-qfMiQ,-q)-' 

where Q = e ~* with t the Kahler modulus of the base P^ "^-t- X such that \Q\ < 1, and 

is the generalized MacMahon function, the generating function for weighted plane partitions. 
The extra parameter Q weights the contributions from "fractional instantons" stuck at the res- 
olution of the conifold singularity, or D2-branes in the D6-brane gauge theory on X, with charge 
equal to the number of diagonal boxes of the plane partition. On the other hand, the Donaldson- 
Thomas theory of the noncommutative crepant resolution of the conifold singularity enumerates 
framed cyclic representations of a quiver with the Klebanov-Witten superpotential [53], or equiv- 
alently of the corresponding quotient path algebra A. The corresponding partition function can 
be computed by counting partitions of a two-coloured pyramid of length one|3 associated with 
a perfect matching of a brane tiling, and is given by [53] 

(2.11) Z^^{q,Q)=M{-qfM{Q,-q)-'M{Q-\-q)~' . 

The appropriate measure in this case involves the Hall-Littlewood polynomials px{x;v) [20] , 
one of two generalizations of the Schur polynomials. They are symmetric polynomials in x, 
homogeneous of degree |A|, and with coefficients in Z[w]. As with the Schur polynomials, the 
Hall-Littlewood polynomials can be extended to Hall-Littlewood functions involving an infinite 
number of variables. They are defined by 

(2.12) p,(.;.)= ^ a{x J] ^ 

aeS„/S>>, A,>Aj ■' 

where S^ is the subgroup of the symmetric group Sn consisting of permutations that leave 
A = (Ai, . . . , A„) invariant, and a (/(x)) = / (o" {x)) with (j{x) = x"^ ■ ■ ■ x^". The parameter v 
serves to interpolate between the Schur polynomials 5\{x) at f = and the monomial symmetric 



Other 5-deformations are introduced in [14] . where only one of the Schur functions is g-specialized. 
The two colours are weighted by the variables go = g/Q and gi = Q. 
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functions mx{x) = J2a£Sn/s^ '^(^) at v = 1. The corresponding Cauchy identity for Hall- 
Littlewood polynomials is 

(2.13) J]6A(t;)pA(x;t;)pA(2/;t;)= J] ^r^^^^^ = ^hl , 
where 

Ai 

(2.14) bxiv) = llivU^(x)- 

i=l 

Here {v)o = 1, {v)m = ni<j<m (^ " ^■') ^^ m G N, and {v)oo = HjeN (l " ^■')' "^^ile mj(A) is 
the number of parts of A equal to i. 

In parallel to the C^ case above, if we make the specification Xi = (— g)*~^'^ and yj = (— ^)-'~^'^ 
in ()2.13p . and set the parameter v equal to either Q or Q^^, then we obtain expressions for the 
conifold Donaldson-Thomas partition functions analogous to (12 .Sp given by 

(2.15) Z$T {q,Q) = M{-q) ^ hx (Q) Pa( {-qf-^'^ ^Qf 

A 

and 

(2.16) Z^^{q,Q) = Y,bx{Q) Px{{-qr''^ ;Qf ^ ^ (Q"') P,{i-<ir"^ '^Q" 



A M 

Likewise, one can define a Hall-Littlewood measure on partitions given by 



(2.17) A^HL {A} = ^5^ hx{v) Pa(x; t;) pA(y; t;) , 

with the normalization 2^hl given by (j2.13p . This leads to the Donaldson-Thomas partition 
functions on the conifold ()2.15p and (I2.16p . when specializing to Xj = (— g)*~^" and yj = 
{—qy~'^''^- This specialized measure can be regarded as a Hall-Littlewood deformation of the 
(7-Plancherel measure. 

We have seen that the generalization of symmetric functions provided by the Hall-Littlewood 
polynomials is appropriate to the Donaldson-Thomas partition functions of the conifold, in con- 
trast to the simpler case of C^. This is especially true in the noncommutative setting (|2.16p . 
where even the two MacMahon function factors in (|2.1ip are automatically included. The ap- 
pearance of Hall-Littlewood polynomials in topological string theory has already been discussed 
in [24]. In particular, the expression (|2.15p for the partition function of the resolved conifold is 
given there. However, the reduced partition function 

(2.18) Z^r (q, Q) = %^ = MiQ, -q)-' 

factors out the contributions from the degree zero subschemes (regular DO-branes) of X, and it 
can be written solely in terms of Schur functions by using the dual Cauchy identity 

(2.19) ^SA(x)sA'(y)= n i'^ + ^iVj) 

together with the scaling property 5x{Qx) = Q' ' 5x{x) to write 

(2.20) Zd\ (g, Q) = Y: (-Q)I'I Sa ii-qy-'^') ^A' {{-qy-'^') • 
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Similarly, one has 

(2.21) 

X,fi 

2.3. Orbifold. Consider the action of the cyclic group Z2 on C'^ generated by {zi^Z2,z^) h^ 
(— zi, — 2:2, -^3)- The crepant resolution of the orbifold singularity C^/Z2 given by the Z2- 
Hilbert scheme is Y = C'pi(— 2) © C'pi(O) = Opi(— 2) x C, where the first factor is the min- 
imal (Hirzebruch-Jung) resolution of the Ai Klein singularity C^/Z2 [21]. The corresponding 
Donaldson-Thomas partition function is similar to that of the conifold ()2.9p and reads [54] 

(2.22) Zl^{q, Q) = M{-qf M{Q, -q) , 

where again Q = e~*. The Donaldson-Thomas invariants of the noncommutative crepant 
resolution, given by the quiver algebra of the McKay quiver associated to the Ai singularity |21] . 
are the same as the orbifold Donaldson-Thomas invariants of the quotient stack [C^/Z2] [Ml 
155] . They are labelled by Z2-representations p and enumerate Z2-invariant zero-dimensional 
subschemes y C C^ with H^{Oy) = P- The corresponding partition function can be computed 
by counting configurations of two-coloured boxes, with the colours corresponding to the two 
irreducible representations of Z2 , and is given by 



(2.23) Z^f\q,Q) = M{-qf M{Q,-q)M{Q-\-q) . 

The reduced orbifold partition functions are "dual" to those of the conifold, in the sense that 
they are related through 

(2.24) Z^r,{q, Q) = Z^^{q, Q)"! and Z^t^^^(q, Q) = Z^Q, Q)"' • 

Hence one can regard the orbifold generating functions in terms of Hall-Littlewood measures by 
using the formulas of the preceding subsection. Alternatively, they can also be written solely in 
terms of Schur functions as 

(2.25) z^T{q,Q) = E Q'"sa((-?)^-^/')' 

A 
and 

(2.26) Zj?^/"^ {q,Q) = Y: Ql^l-l^'lsA((-g)^-^/^)^.((-g)^-^/^)' . 

X,fi 

2.4. BPS state counting. For the toric Calabi-Yau backgrounds X of interest in this paper, 
which have no compact divisors, it was shown in [561 [57] that the closed topological string 
partition function can be expressed in terms of genus zero Gopakumar-Vafa invariants as 

(2.27) Z4(g, Q) = M(g)>^(^)/2 J] J] M(g- Q^ g)"^™ 

l3eH2{X,Z)+ m 

where x{^) is the topological Euler characteristic of X, the integer Nm counts the number 
of BPS states of M2-branes in curve class /3 = (/3i, . . . ,/3s) G H2{X,1i) and with intrinsic 
SU{2) spin m, and Q^ = Q^ ■ ■ ■ Qs . This expression describes the generating functions Z^pg, 
counting BPS bound states of a single D6 brane wrapping X with D2 and DO branes, in different 
chambers of the Kahler moduli space which are separated by walls of marginal stability where 
BPS states decay or form. This includes the chambers that lead to the Donaldson-Thomas 
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partition functions at both commutative and noncommutative points in the Kahler moduli 
space [531 EH [59] . Iii any chamber, the BPS partition function is given by 

(2.28) Z^PS = ^4(q, Q) Z,%iq, Q-i) . 

The symmetric function expansion in this case involves the other generalization of the Schur 
polynomials, the Jack polynomials j\° (x) [601 [20] with j^"~ (x) = S\{x). They satisfy a Cauchy 
identity 

(2.29) E^A(a)ji"^(x)jl")(?/)= n (l-^.y.)"'/" , 

where cx (a) are rational functions of the parameter a which have been calculated in [60]. Due to 
the power 1/q in ()2.29p . they are well suited to describe the topological string partition function 
as written in ()2.27p . The generic expression (|2.27p can thus be written in terms of the Jack 
polynomials as 

^'■''^ w§(Sj^ = n n E ^a( - i/Ni) ,-1^1 Qi^i/^ii-v-^)(,-v2)2 

^^' l3eH2{X,Z)+ rn A 

and the powers of the MacMahon function can be expressed as 

(2.31) M(g)x(^)/2 = Y: c,{2/x{X)) j(2/xW) (,-1/2)2 . 

A 

For the backgrounds considered in this paper, the (unreduced) noncommutative partition 
functions are related to those of X by the wall-crossing factor W^ = Z^r^{q,Q^^), which 
describes the crossing of an infinite number of walls in going from the noncommutative point 
to the large volume point [53t [58] . We can index the walls crossed by i G N and factorize the 
wall-crossing factor as W^ = Y[i>i ^j^j with W^^ the jump of the BPS state partition function 
across the i-th wall. This yields the partition function ()2.28p in the n-th chamber 

(2.32) Z]^Ps(n) = W^Z^T(9,Q), 

where W^ = Ui>n W^^ ■ 

For the resolved conifold X = X, this prescription yields W^ = (1 — (—(7)*)* and 

(2.33) w^^ ((_,)" Q-i; _,);^Zd\(9, i-qrQ-') 

where (a;g)oo = JHigN (^ ~ '^Q'^)- The expansion of (j2.32p in Hall-Littlewood functions is given 
by 

A 

(2.34) X J2 K H-qT Q-') p^H-qT''^' ; {-qT Q"')' • 

/^ 

This function interpolates between Zgpg(l) = Z^rj, (after reparametrization) and Z^pa^{oo) = 
Zprp, at the noncommutative and large volume limit points, through values of n corresponding 
to chambers in the Kahler cone of X. For the resolution of the Ai quotient singularity X = Y, 
using duality one has Wj = (W^) . The expansion of the BPS state partition function (|2.32p 
in Schur functions is then given by 

^^ps(^) = {{-qr Q-' ; -q)1 

(2.35) X Y: Q'^''"'^' i-qr^''^ Bx{{-qy-'/')' B,{i-qy-'/'f B,{{-qy-'/')' . 
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2.5. Generalizations. These constructions generalize in a natural way to all local toric Calabi- 
Yau threefolds without compact divisors. Let us briefly illustrate this through a representative 
set of examples, generalizing the conifolds and Z2-orbifolds considered above. Let G be a fi- 
nite subgroup of SU{2) acting on C^ via the natural embedding SU{2) C SU{3). The crepant 
resolution of the orbifold C^/G given by the G-Hilbert scheme of C^ is Yq = Sg x C, where 
§G is the minimal ADE resolution of the double point singularity C^/G. By the McKay cor- 
respondence |2Ij, the non-trivial irreducible G-representations correspond to simple roots of an 
associated ADE root system A, the collection of which gives a basis for H2{Yg, Z)- Let A+ C A 
be the set of positive roots, and n the number of irreducible representations of G. Then the 
Donaldson-Thomas partition function generalizes ()2.22p to ^61j 

(2.36) Zl^{q,Q)=M{-qr f] ^^^ -<?) ■ 

/3GA+ 

The noncommutative Donaldson-Thomas invariants arising from the McKay quiver associated 
to the afiine ADE Dynkin diagram are the same as the orbifold Donaldson-Thomas invariants of 
C^/G [55l|6T]. The corresponding partition function Z^^r^ {q,Q) is given by the same formula 
()2.36p but with the product now ranging over the full root lattice A |61] . Hence the corresponding 
reduced partition functions can be expanded in Schur functions as 

(2.37) Z^^{<1,Q)= n E Q""^A((-g)^-^/^)\ 

/3eA+ A 

and the identical formula for Z-^rj! (q, Q) involving a product over all roots. The corresponding 
BPS state partition functions in the various chambers can be worked out as before. 

For the special case G = Z^, the products run over roots of the An-i Lie algebra and the 
partition function (|2.36p specializes to [Ml [53] 

(2.38) zl^^-{q,Q)=M{-qr J] M{Q[,^^^,-q) , 

l<i<j<n 

where Q[ij] = Qi Qi+i ■ ■ ■ Qj- This is (after reparametrization) the generating function for n- 
coloured plane partitions. It is "dual" to the Donaldson-Thomas partition function for the 
generalized conifold geometry X„ \Q2\ [63| [56] which is given by 

(2.39) Zl^{q,Q) = M{-qr \{ M(Q[,,,], -(?)^- , 

l<j<j<n 

where Nij = — (— 1)"»j with riij the number of internal edges between vertices i and j in the 
toric web diagram for X,^. The corresponding reduced function is expanded in Schur functions 
as 



Z^^iq.Q) = n E(^M)'"sA((-g)'=-^/')' 



l<i<j<n A 

Nij=+1 

(2.40) X n E(-Q[M])''^'v((-^)'^'/')v((-9)'"'/')- 

l<j<j<n A« 

These expressions should all follow from the representions of the Donaldson-Thomas partition 
functions as correlators of vertex operators given in |64) . as it is known that such correlation 
functions can be represented in terms of Schur functions. These formulas also formally apply to 
the C^/Z2 X Z2 orbifold and its symmetric resolution, the closed topological vertex geometry. 
The noncommutative chamber was originally considered in |54) . while several infinite families of 
chambers and their wall-crossing behaviours are identified in [Bl]. They can also be extended 
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to some non-toric geometries without compact divisors, such as the formal toric Calabi-Yau 
threefolds considered in |65t [66] (see also [56] ) . The resulting expressions are straightforward but 
somewhat tedious to write down, and one must be somewhat more careful with the convergence 
of the various infinite products involved. We omit the details. 

3. Unitary matrix models for Donaldson-Thomas theory 

3.1. Chern-Simons theory and Toeplitz determinants. The interest of Toeplitz determi- 
nants in physics begins with Onsager's study of the two-dimensional Ising model, as he succeeded 
in showing that the diagonal spin correlation function is given by an A^ x A^ Toeplitz determi- 
nant |67[ [68] . Another classic example is the one-plaquette model of two-dimensional Yang-Mills 
theory with gauge group U{N), whose partition function is given by |69j 

(3.1) Z^(A)= det [/,„,(2A)] , 

where Iniz) is the modified Bessel function order n. The same problem was shortly afterwards 
studied using unitary one-matrix models jTD], leading to the well-known Gross- Witten matrix 
model 

^ /■27r j/i . , 

(3.2) Z^(A) = n / 5exp(A(e^^^ + e-^0) n 

j=i -^o k<i 

In recent work [31] we showed that the unitary matrix model for U{N) Chern-Simons gauge 
theory on the three-sphere S^ is, in a certain sense, a g-deformation of the Gross- Witten model. 
Here we shall show that its partition function can also be represented as a Toeplitz determinant. 
The reason for the existence of two equivalent representations is a direct relationship between 
Toeplitz determinants and unitary random matrix models, found long before the introduction 
of matrix models. It is given by the Heine-Szego identity [26] 

N I-2TT J/1 

(3.3) n / ^fi-''') nh''^-«^'f =^^(/)' 

j=i -^o k<i 

where D]\j{f) is the N x N Toeplitz determinant with symbol /, 

(3.4) Z);v(/) = ^ det \f{j-k) 

l<j,k<N L 

with /(r) = 2^ fn^ f{^^^ ) d^, r G Z the Fourier coefficients of the symbol function. The 
symbol function f(z) is the weight function of the corresponding matrix model; we assume 
throughout that it lives in C°°(S'^). In the case of the Gross- Witten model, one uses the identity 

oo 

(3.5) exp(A(z + z-i)) = J^ /„(2A)z" 

n=— oo 

to establish the equivalence of the two expressions ()3.ip and ()3.2p . 

Let us now compute the U{N) Chern-Simons free energy on S^ in terms of a Toeplitz 
determinant. The unitary matrix model is given by the partition function |29] 



(3.6) ^cf^(5^)-n/^^0(«^'^i«)nh^' 



e 



j = l •'0 ■• fe<; 



with 



(3.7) e{z\q)= Yl ^' ^'^' 
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We have seen that the entries of the Toephtz determinant (|3.4|) are given by the Fourier coef- 
ficients of the symbol, which is the weight function of the unitary matrix model. In this case, 
the symbol is the theta-function Q{z\q) and the entries of the matrix are automatically at our 
disposal. Then the determinant is 



(3.8) DNiO) 



ao ai • • • QN-i 

0.-1 ao ■ ■ ■ aN-2 

O-JV+l a-N+2 ■ ■ ■ oo 

2 , 



■2 /o 

with Un = a-' ' 



Hence the required determinant is detj j [q (* ^' ' ^] , which is known to give rise to the Chern- 
Simons partition function [40] 

Af-l 

(3.9) Z^TH^^) = ,iet, \^-^-^^^'] = n (1 - 'i^-' ■ 

- '^- k=i 

Notice that the non-trivial part of the determinant, detj^- [(?*■'], is just a Vandermonde determi- 
nant. 

A more general Toeplitz determinant, with symbol 



(3.10) Yl 



oo 

2 



dV2 J 



Z' 



j=0 

was already computed in [71j. It was proven there, by induction, that the determinant ()3.8p 
with coefficient am instead of oq is given by 

N-l 

(3.11) Dm,N = {q-^^^-^yy n (1 - 'iT~' ■ 

k=l 

They consider the case m > N because their symbol is (|3.10p . but if the symbol is the full 
theta-function ()3.7p we can take m = and recover Dq^n = Z^^ (S^) given in (|3.9p . 

Of course, knowledge of the Toeplitz determinant result directly implies, due to p.3p . that 
there is a unitary matrix model representation with a theta-function as weight function, a model 
found by other means in [29]. In [30] it was shown that the inverse of the theta-function also 
leads to a viable Chern-Simons matrix model. We shall see that there is an equivalence between 
these results and the ensuing representation of Donaldson-Thomas partition functions in terms 
of unitary matrix models. 

3.2. Donaldson-Thomas theory on C^. As in the case of Chern-Simons gauge theory on 
5^, we can write the partition function of Donaldson-Thomas theory on C'^ as a Toeplitz deter- 
minant. For this, we use Gessel's formula for the product of Schur polynomials in terms of a 
Toeplitz determinant, which in terms of the Schur measure reads [25] 

(3.12) VNix, y)= Yl ^^ (^) ^^ (y) = Dn{A-j) , 

X:Xi<N 

where 



(3.13) Ai = Ai{x, y) = y] i)i+iix) Uy) 



oo 

E 

1=0 



and f)r(a^) = X^i ^I i^ ^^^ r-th complete symmetric function. Then the Donaldson-Thomas 
partition function is given by 

(3.14) zSt (q) = lim Pjv(x. = q'-'^' , y. = q'-'/^) . 

Af— 5>oo 
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It is convenient to consider the symbol of the Toephtz determinant [25\ [72] 

oo 

(3.15) /c3(^)= Yl Mx,y)z' = lli'^-y^^'T' (^-^^^y' ■ 

i=—oo j>l 

Taking into account the Heine-Szego identity (j3.3p . the Donaldson-Thomas partition function 
can then be written as the A^ — )• oo limit of an A^ x A^ matrix model 



n=l -^^ j=l 



(3.16) X JJ le'^*" - e 



^ieH2 

k<l 

Using the Jacobi triple product formula for the theta-function 



(3.17) Q{z\q) = 11 {I- q') (l + q'''^^ z'') [l + q^'^'^ 



we can write 

(3.8, .c.(,^,^nf^ep^n 



e '' — e 



where Aqo = limTv-^-oo (9)00- Thus the Donaldson-Thomas partition function is a unitary N = 00 
one-matrix model with weight function w {9) = 0(— e ^^\q)~^- 

Recall that the matrix model for U{N) Chern-Simons gauge theory on S^ is a unitary matrix 
model with weight function w' (9) = 0(e '^|g) and N eigenvalues |29[ [3T] . We relate the two 
models and the corresponding partition functions in a more precise way below. Before doing that, 
we can also demonstrate the relationship using an equivalent matrix model representation of the 
Donaldson-Thomas partition function, based again on the Cauchy identity but now written as 

(3.19) Yl ^A' (x) sv {y) = n T^^— 

where A' are the conjugate (transposed) partitions. The symbol in this case is [^ 

(3.20) /^3 (z) = n (1 + yj ^"') (1 + ^^- ^^ 

and the corresponding matrix model representation is 



Uj ^1 , 



/"ZTT J/1 



j-1/2 ^u 



n=l -^0 ■■ i=l 



(3.21) X n 



e^^^-e^^f 



Using ()3.17p again, the expression in terms of theta-functions is then 

(3.22) <(<?) = Coo n r ^ Qi^''-\l) WW" - ^'''W 

n=l -^0 ^ ^^; 



-Af 



with Coo = liniAf-i-oo (Q')oo • The matrix model representation of the MacMahon function [50] 
has not appeared before. 
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3.3. Chern-Simons/Donaldson-Thomas correspondence. The two matrix model repre- 
sentations, ()3.18p and ()3.22p . parallels the situation in Chern-Simons theory, since the weight 
function w {9) = Q{—e^^\q)~^ can also be used for a Chern-Simons matrix model. The Chern- 
Simons unitary matrix model representation thus not only matches (|3.22|) but also (|3.18|) . This 
result appeared in j30[ eq. (4.8)], where the non-uniqueness of the matrix model representation 
was pointed out, based on the undetermined moment problem for the Stieltjes-Wigert weight 
function. It follows from the explicit expression, given by Askey, for a weight function that 
has the same moments as the log-normal distribution, which is given by the inverse of a theta- 
function [73] 



If we set z = y^ e '^ and choose 7 = 0, then by the triple product formula ()3.17p we directly 
obtain the Chern-Simons — or finite N — matrix model version of ()3.18p n 

Thus in both cases, we see that the right-hand side is essentially the matrix model for U{N) 
Chern-Simons gauge theory on S"^ when N = oo. Indeed, considering the corresponding gener- 
ating functions, we have 

(3,24) Z,T'(S')=n^^ a„d Z&(,) = n^. 

Then I^T22^ implies Zg^ (q) = Coo liniAr^oo Z^^^\s^), which also follows from I^Tm and 

oo 

(3.25) iirJ = 11(^-1')''' ■ 



3 



Prom the very definitions of their partition functions, the Donaldson-Thomas free energy F^rj, 
log Zprp and the Chern-Simons free energy Fq^ = log Z^c, at A^ = oo thus satisfy the simj 
relationship 



oo 



(3.26) F§'^ (q) = F^t\s') - lim N ^ log (l - q^) . 

i=i 
Hence, via a precise infinite renormalization, the Donaldson-Thomas free energy is exactly the 
free energy of U{N) Chern-Simons gauge theory in the limit N — ?■ oo. 

The factor liniTv-s-oo (q)^ that links the two partition functions has a natural interpretation. 
As we have seen, the g-Pochammer symbol 

n 

(3.27) (g)„ = n (1 - ^') 

is the building block of the Chern-Simons partition function. In the limit n — )• oo, it also appears 
as the partition function of U{N) topologically twisted Vafa-Witten A/" = 4 gauge theory [71] 
on C^, which is given by [21] 

(3.28) ZC;^) (q) = (g)-^ . 
We can thus write 

(3.29) zSt (q) = lim Z^J^v) (q) Z^Hs') . 

A/— ^oo ^ ^ 



In [3D], the particular case 7 — —3/2 was considered, which gives one of the most common normalizations of 
the Stieltjes-Wigert polynomials. 
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It would be interesting to better understand this relationship from the point of view of the six- 
dimensional U{1) topological gauge theory underlying the Donaldson-Thomas invariants OH]. 
If we regard C^ as the trivial line bundle over C^ , then this result suggests that the fibre degrees 
of freedom in the six-dimensional gauge theory can be integrated out, leaving its natural gauge 
theory counterpart on C^ and on the boundary S^, at infinite rank. This sort of reduction of 
the gauge theory partition function is demonstrated in the case of local toric surfaces in [75] . 

However, in (|3.28|) the (complexified) gauge coupling parameter is q = exp(27rir) with r = 
-^ (for vanishing 9-angle). If we use the relationship with the topological string coupling 

9ym 

constant gs = g\y[/2, then the quantum parameter is q = exp(— 47r^/(7g), while the quantum 
parameter in Chern-Simons theory is q = exp {—gs) ■ Thus in ()3.29p the Chern-Simons partition 
function should be understood in its dual form after performing a Gauss summation. This is 
the correct setting for merging the four-dimensional Vafa-Witten theory with its corresponding 
boundary Chern-Simons gauge theory [76l[75]. Indeed, the equality ()3.29p can be regarded as 
the large A^ limit of the partition function for g-deformed Yang-Mills theory on the two-sphere 
5^ [76J, which describes the relationship between the two gauge theories. We will discuss the 
relation between Donaldson-Thomas theory and g-deformed gauge theories in the next section. 

3.4. Conifold Donaldson-Thomas theory. The case of the conifold partition functions can 
be similarly studied. Using sx{Qx) = Q''^'Sx{x), the symbol associated to the slightly more 
general expansion (|2.20p is given by 

(3.30) fx{z) = ll{l-QyjZ-yU^ + Xjz) ■ 

i>i 

Hence the Donaldson-Thomas partition function of the resolved conifold (without degree zero 
contributions) has the matrix model representation 

(3.31, zs^M^-^^^Yii lfn .i;;-.;:-.».. nh'"-e'"r 

n=l " j=l ^ ^ k<l 

The noncommutative conifold partition function (|2.2ip can be written as the product of two 
Toeplitz determinants !?(/) = detr(/), where r(/) = [/(j — A;)] . ^^^ is the corresponding £^(N) 
Toeplitz operator. However, the product of two Toeplitz operators is not of Toeplitz form |77j . 
To deal with this problem, one considers all Toeplitz matrices T/v(/) = H^v T{f) Hjy for N G'N 
together, where H^r is the orthogonal projection onto the Fourier modes 1,...,A^. Then a 
sequence of products of Toeplitz matrices {Ttv (/) T/v ig)}^^!^ ^^ asymptotically equivalent to 
the sequence of Toeplitz matrices {T^r (/5')}jveN t^- "^^^ requisite condition is that the Fourier 
coefficients of the symbol tk = f{k) (the entries of the Toeplitz matrix) are absolutely summable 

(3.32) Yl \*k\ < oo . 

fc=— oo 

A sequence of Toeplitz matrices T^r = [tj-k]i<j,k<N for which the tk are absolutely summable 
is said to be in the Wiener class. Likewise, a function / (z) defined on the circle S^ is said to 
be in the Wiener class if it has a Fourier series expansion with absolutely summable Fourier 
coefficients. If the symbols are of Wiener class, then one has the convergence result [TTj 

(3.33) hm ||r^(/)r^(r7)-r^(/5)||=o 

where || — || denotes the operator norm on finite-dimensional matrices. 

In our case, the symbols are given by theta-functions which have Fourier coefficients of the 
type tk = q^ '^ with < q < I. Absolute summability is consequently satisfied and hence the 
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symbols are of Wiener class. The limit A^ — >■ cx) is precisely the one that we are studying, and 
hence we can write 



^ .2.^00 (1 + ^.-1/2 eie„)2 

27r J--I- 

n=l "^ j=l 






(3.34) X JJ I e ''- - e 



^ 9fc _ ^ie,|2 
with the limit understood in the sense of norm convergence. 



3.5. Orbifold Donaldson-Thomas theory. The orbifold partition functions are treated anal- 
ogously. The symbol associated to the Schur function expansion ()2.25p is 

(3.35) /y (z) = J] (1 - Q yj z"^) "' (1 - x, z)-^ . 

This leads to the matrix model representation 

N ,2-K 



ZD\(,,g) = lim n / ir n n n . 1 



/2 e-i^") (1 -gi-1/2 gi^-) 



(3.36) ^ n I ^ ' 



?fc _ e'^'l^ 



fc<« 



Likewise, by using the same arguments which led to ()3.34p . the noncommutative orbifold parti- 
tion function ()2.26p admits the matrix model representation 

^ '■2-d^- (l_g:;-i/2 ei^")-" 

.„ 27r J-i 

n=l " " i=\ 

(3.37) X n 

k<l 



Z'^''^^ (a O) - lim TT / ^ n (1-g^'^^^^^") 

DT W.t^J ^irn^j^l/^ 27r y (l-Qgi-1/2 e-i^") (l-Q-igi-V2 e-ie«) 



e^^^-e^^f 



As pointed out in [7^, the arguments which led to (|3.34p easily extend to the more general 
case of products of m, > 2 Toeplitz matrices. In particular, one has the convergence results 

(3.38) lim \\Tn (/i) ■■■Tn {fm) - T^ (/i • • • fm) || = . 

This enables us to write down matrix model representations for all generalizations considered in 
the previous section. For example, the BPS state partition function (|2.35p can be expressed as 

^ r27r JO 



4psW = {{-qrQ-';-q)lC^ V^m f[ [ 



2tt 

N-2 



f3 39) X TT [l-q' I e , TT I e '^^ - e '^'l' 

^ ' l^i (1-Q (7^-1/2 e-i^-) (l-Q-i(-g)"gi-V2 e-i^-) y ' 

Similarly, the generalized orbifold partition function (j2.37p is 

(3.40) 

N .27r J/1 ~ (-, _„j-i/2 pie„^-l 

ZJ«(g,Q)= lim TT / — n TT ^ ^, . J, ^ TTh'^^-e'^f , 
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while the generahzed conifold partition function (|2.39p can be expressed as 






m=l -"^ p=l l<i<j<n "*l«JJ^ 



e'^f 



and so on. 

The Toephtz determinant can also be expressed as a Fredholm determinant [78} [79] , which 
suggests a deep underlying integrability structure, as we discuss in detail in the next section. 
This applies directly to the Chern-Simons partition function, which is given by a finite N x N 
Toeplitz determinant, for which the result of |79j immediately applies. This result does not 
strictly apply in the limit A^ — )• oo which yields the Donaldson-Thomas partition function. In 
this limit, the Fredholm determinant representation converges to 1 and the infinite Toeplitz 
determinant is given by the normalization constant of the Schur measure [79]. This agrees 
with our computation, wherein the matrix model representation of Donaldson-Thomas theory 
follows directly from the Toeplitz determinant representation of the Cauchy identity (|3.12p and 
the Heine-Szego identity ()3.3p . Proper specialization of the symmetric functions then yields the 
Donaldson-Thomas partition functions. For the first part of this procedure, one can alternatively 
use the results of [80], where the Cauchy identity is expressed directly as a generic matrix model 
average. Other instances where the partition function can be written as a Fredholm determinant 
are the partition function of two-dimensional quantum gravity [81] and the grand canonical 
partition function of c = 1 string theory with vortex excitations [82] . 

4. Integrability structure 

4.1. Toda and Toeplitz lattice hierarchies. The expansions in Schur functions are also use- 
ful to establish a relationship between Donaldson-Thomas theory and the theory of integrable 
hierarchies, a connection that is generally expected whenever a matrix model formulation is avail- 
able |83[ ^^ . Using the results of [85] , it is straightforward to identify the Donaldson-Thomas 
partition functions as particular instances of a tau- function of the 2-Toda lattice hierarchy |86] . 
More precisely, they are given by tau- functions of the Toeplitz lattice hierarchy |85[ 187] . which 
is a reduction of the 2-Toda lattice hierarchyjj The tau-functions of the 2- Toda lattice hierarchy 
Tn{t,s), n G Z depend on two sets of time variables t,s £ C°° and are defined by the Hirota 
bilinear equations. They can also be written as the determinant of a semi-infinite moment matrix 

(4.1) Tn{t, s) = det TJlnit, s) . 

In the case of the reduction to the Toeplitz lattice hierarchy, this moment matrix is a Toeplitz 
matrix (whereas in the reduction to the standard Toda lattice hierarchy it is a Hankel matrix), 
and the components of the vector T{t,s) = (To(t, s) = l,ri(t, s), . . . ) of tau-functions of the 
Toeplitz lattice hierarchy satisfy [85] 

(4.2) Tn{t,s)= Y, 5x{t)5x{-s) . 

X : Ai<n 



As pointed out in [85], the Toeplitz lattice hierarchy is better known as the Ablowitz-Ladik hierarchy [88] 
which arises in discretizations of the non-linear Schrodinger equation, but we follow the terminology of |85l 187] . 
The other reduction of the 2-Toda lattice hierarchy leads to the standard Toda lattice hierarchy [SS] . The Toeplitz 
reduction (on S^) and the Toda reduction (on R) are essentially equivalent |89| . 
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Hence if the times are taken to be ti = g*~^/^ and Sj = —q-'~^''^, then 

(4.3) r^ir'^\-q^-'/') = zE'^{q) . 

On the other hand, if the two sets of time variables are taken to be ti = ^/Q q^"^''^ and Sj 
-\/Qg^'"^/^ then 



(4.4) r^iVQq'-'/^-VQq'-'^'') = Z^^{q,Q) = Z^^{q,Q)-\ 

Thus while one begins within the formalism of the 2-Toda lattice hierarchy, the Donaldson- 
Thomas partition functions are a particular case of a reduction of the 2-Toda lattice hierarchy. 
This reduction is equivalent to the one-dimensional Toda lattice hierarchy, a result that has also 
been shown in the context of the melting crystal picture [23]. While our partition functions 
involve two sets of Schur polynomials, the two sets of times are taken to be equal. That this 
reduces the system to a one-dimensional Toda hierarchy is shown in [23] . Using the free fermionic 
representation of a deformation of the expansion ()2.25p . they show that it is a 2-Toda tau- 
function that satisfies 

(4.5) r(t,s) = r(i-s,0) , 

and thus reduces to a one-dimensional Toda lattice hierarchy. 

More generally, the tau-function of the two-component KP hierarchy is also a function of two 
sets of time variables and admits the double Schur function expansion [32] 

(4.6) rKp(t,s) = ^ cx^,5x{t)5^{-s) , 

A, /I 

where c^^ are Pliicker coordinates of a point on a two-component analog of the Sato grassmann- 
nian [32J and have well-known determinant expressions. (The Schur expansion for tau-functions 
of the ordinary KP hierarchy consists of only one set of Schur functions.) The case of the 2-Toda 
lattice hierarchy follows from this one, as it only has an additional discrete index n as in (j4.2p . 
More precisely, following [90j these tau-functions can be written as images of the Pliicker map 
corresponding to projection along a basis element of the charge N sector of a fermionic Fock 
space as 

(4.7) r^Z^t.s) = Y, BM,g[\tA sxit)sx{-s) . 

\,fi 

Here Bjy^g (A,/x) are Pliicker coordinates of the image of an element g G Gr^N{H) in the grass- 
mannian of subspaces of the Hilbert space % = L?'{S^) whose orthogonal projections onto the 
subspace ?^+ C %, consisting of functions that admit holomorphic extensions to the interior 
of S^ C C, have Fredholm index N . Then the noncommutative Donaldson-Thomas partition 
function Z-^r^ {q, Q) is a particular case of (|4.7p with unit Pliicker coordinates. 
The particular case of (|4.7p for the identity element g = e\s [90] 

(4.8) 4^i(t,.) = j;sA(t)sA(-.) 

A 

which, with the specification of the time variables given above, describes the partition functions 
■^DT (9) ^'^'i Z^r^{q,Q), and reduces to a one-dimensional Toda lattice tau-function. With 
the matrix model expression for Z-^r^{q,Q), one can interpret it as the product of two one- 
dimensional Toda tau-functions. The more general partition functions described earlier are 
likewise given as products of 1-Toda tau-functions. 
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4.2. Isomonodromic tau-functions. Following [91J, a somewhat more conjectural relation- 
ship with tau-functions and integrable systems can also be established, using the fact that the 
Donaldson-Thomas and Chern-Simons partition functions are Toeplitz determinants. The result 
of [91j shows that some Toeplitz determinants can be interpreted as certain isomonodronic tau- 
functions, which coincide with those of Jimbo, Miwa and Ueno [92j . They focus on a particular 
case of the Gessel identity (|3.12p where one set of variables x = {xi, . . . , x„) is generic and left 
unspecified, while the other set has the specification y = (1, . . . , 1). In this case the associated 
symbol is 

n 

(4.9) hMv{z)=e^''\{{l + x,z) . 

i=i 

However, as pointed out in [91], their derivation of integrable partial differential equations can 
be applied to any Toeplitz determinant for which the logarithmic derivative of the associated 
symbol f{z) is a rational function of z. 

The Chern-Simons partition functions for gauge group U{N) have finite-dimensional Toeplitz 
determinant representation and symbol given by a theta-function ()3.17p after the proper spec- 
ification, i.e. n — 7- oo. This symbol violates the rational behaviour required of its logarithmic 
derivative, unless we truncate the products at some large value M < n. The extension to infinite 
support of the variables x = (xj) is left as an open problem in [91j. In this "finite-dimensional" 
approximation, the calculation of the Toeplitz determinant in Section [3.11 proceeds based on the 
Fourier coefficients up to au-i, and corresponds to a truncation of the symbol to index M. 
Then the truncated Chern-Simons partition functions can be identified with the Jimbo-Miwa- 
Ueno tau-functions corresponding to the (generalized) Schlesinger isomonodromy deformation 
equations of the 2x2 linear matrix ordinary differential equations which have M simple poles 
in C and one irregular singular point at infinity of Poincare index one. 

The rigorous proof that the Chern-Simons partition function is an isomonodromic tau-function 
is beyond the scope of the present paper L_| It is then an interesting open problem to show that the 
Donaldson-Thomas partition functions are tau-functions of a similar system of non-linear partial 
differential equations, which have both infinitely many poles in the complex plane and are given 
by infinite-dimensional Toeplitz determinants. The isomonodronic tau-functions are known L93j 
to be intimately related to the KP tau-functions discussed above and introduced in [32], hence 
it should be possible to further establish a stronger relationship between these two fainilies of 
integrable hierarchies and Donaldson-Thomas theory. The expressions of generic Donaldson- 
Thomas partition functions as correlators involving exponentials of fermionic bilinears given 
in [63] shows that they automatically provide tau-functions of integrable hierarchies. 



4.3. Nekrasov functions. The expansion of Donaldson-Thomas partition functions in terms of 
the g-Plancherel measure is also useful for establishing a relationship with J\f = 2 supersymmetric 
gauge theory in four dimensions with Casimir operators, since Nekrasov's formulas |94j can be 
also written in terms of the Plancherel measure |im [TT] . This sort of relationship has already 
been noticed, within the crystal melting picture in |22^ [23] and also in [95j within a slightly 
different context. In particular, for the noncommutative U{1) gauge theorM I the instanton 



This proof presumably follows the arguments indicated in ,28, that the essential characteristics are indepen- 
dent of the cutoff M on the length of the representation associated to A. 

In general, one considers an abelian Af — 2 gauge theory with instantons. The noncommutative gauge theory 
is one possibility, whereby the field theory is embedded in the il-backgound with toric deformation parameters 
(ei, £2). Other possibilities are the gauge theories of fractional D3-branes at an ADE singularity or of the D5/NS5 
brane system wrapping a P^ in K3 [111 . 
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partition function Zinst{0',t,ei,e2) can be expressed as a sum over partitions [TUl [H]- In the 
Calabi-Yau case when ei = —62 = h, one has 

fA tn\ 7 f + t^\ X^ (diniA)^ / 1 >r^ , chfc+i (a, A) 

(4.10) Zinst(a, t,h) = \^ (_^2)|A| «^P ( - ^ 2^ *fc k + l 

where t = {tk) are couphng constants, ch^+i (a, A) are the Chern characters of the partition 
A, and a is the vacuum expectation value of the vector multiplets. Thus it coincides with the 
principal specificatioro of the Schur functions 

(4.11) dimA=lim lim 5x{l,q, ■ ■ ■ ,q^~^) ■ 

As in [ini Hi], let us consider the case ii / and ^2 = ^3 = • • • = 0, i.e. all higher Casimir 
operators are turned off. Using the Chern character ch2 (a. A) = a? + 211? |A|, in that case we 
have 

(4.12) Zinst(a, h,n) = Y^ ^"|^^^2)|Ai ^' ^^P ( - ^ "' - ^1 1^1 

A 

Consider now the C^/Z2 orbifold Donaldson-Thomas partition function (|2.25p . and take its 
g — >■ 1 limit using 

(4.13) limSA(g*-i/2)2 = 57""sA(l,l,...)' 

where q = e"^". Then the partition function Z^r^{—q,Q) in the limit g — )• 1 is the non- 
equivariant limit of Nekrasov's partition function with Q = e"*^ and gs = i^ — )• 0, 

(4.14) lim zM-q = -e-'\Q= e"*^) = e'^^"/^''" Z,^,,{a,t,h) . 

The equivariant case which is directly related with Donaldson-Thomas theory is that of the 
five-dimensional A/" = 1 super symmetric U{1) gauge theory (or K-theory version of the original 
gauge theory), because its instanton expansion involves the q-Plancherel measure instead of 
the Plancherel measure above \12 \ 123) 1 ^ Its partition function coincides with Z^r^{q,Q) and 
Zf^r^{q, Q)"^ given in ()2.25p j22j . As we have seen above and as is also shown explicitly in {23] . 
this partition function is a tau-function of the one-dimensional Toda lattice hierarchy and not 
of the two-dimensional Toda lattice hierarchy. 

4.4. Hurwitz theory. Another gauge theory partition function that can be immediately writ- 
ten in terms of the Plancherel measure is the heat kernel expansion of two-dimensional Yang- 
Mills theory |96[ [97] . On the sphere 5^ = P^ , from this expansion or from the matrix model 
representation [98] the partition function can be written in terms of Schur polynomials as 

(4.15) Z?(f)(Fi)=j;.,(l,...,l)2,^^W, 

A 

where Sa(1, • • • , 1) = dim A is the dimension of the irreducible representation of the U{N) gauge 
group corresponding to A = (Ai, . . . , Aat), and C2(A) is the quadratic Casimir invariant of the 
representation. Hence in the large N limit, whereby the gauge group is U (00) , it coincides with 
the Nekrasov partition function (|4.10p with only the first two Casimir operators turned on, an 
observation already made by [lOJ (see also [95]). 

The large A^ limit picks out the chiral sector of the two-dimensional gauge theory which 
receives contributions from only "small" representations A to the partition function (|4.15p . 



Principal means that it utilizes Schur functions (infinite number of variables), as in Donaldson- Thomas 



theory. 



The string coupling Qs here is then interpreted as the circumference of the compactified fifth dimension. 
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weighted by the hnear Casimir operator. Then the U{oo) chiral partition function is related 
to the Nekrasov function ()4.12p . and one has 

(4.16) Z^i^\FY = lim Z^t {-Q,Q) ■ 

Of course, the q'-deformation of this theory is more intimately related to Donaldson-Thomas 
theory, as then no limit is required. We shall discuss this case in more detail below. We first 
point out that the expansion of the partition function in terms of Schur polynomials allows us 
to apply results of Okounkov |34] and relate the partition functions with tau-functions of the 
Toda lattice hierarchy, while it further enables the study of their role as generating functions 
for ramified coverings of P"^ [,35J . 

Let us begin by briefly summarizing the Toda lattice results of [33]. Let H(i{Ci, . . . ,Cs) 
denote the weighted number of (possibly disconnected) d-fold coverings of P^, ramified over s 
fixed points of P^ with monodromies in fixed conjugacy classes Ci, . . . , C^ of the symmetric group 
'S'dLJ The generating function for degree d coverings of P^ whose ramification over 0, oo G P^ 
can be fixed arbitrarily while the remaining b = s — 2 ramifications are simple is given by 

(4.17) Z^^l{x,y,q,Q)=Y, Q"^-^ E ^a« y. ^d(C^, C,, C(2), • • • , ^2)) , 

d,b=0 ' \a\ = \u\=d " ":>' ' 

times 

where C^ is the conjugacy class corresponding to the partition ^ of d, C(2) is the conjugacy 
class of a transposition, and x^ = Y\^ x^^ for a set of variables x = (xi, X2, • • • )• The generating 
function for connected coverings is the corresponding free energy 

(4.18) F^l (x, y, V, q) = log Z^^ {x, y, v,q) = Y, Q" -^ E ^M ^^ ^I 6 {^^, ^) , 

d,b=0 ' \fj.\ = \u\=d 

where H' ^ (/x, v) denote the double Hurwitz numbers introduced by Okounkov which are the 
weighted numbers of connected degree d coverings of P^ with monodromy around 0, c« G P^ 
given by /i and v, respectively, and b are the numbers of additional simple ramifications. 

The difference between simple Hurwitz numbers and double Hurwitz numbers is as follows. 
Let // be a partition of an integer d. Then the simple Hurwitz number hfi is the number of 
connected genus g branched covers of P^ with monodromy given by fi over a fixed point, usually 
identified with oo, and an appropriate number of fixed simple branch points. In the more general 
setting of double Hurwitz numbers, there is a useful explicit expression for the partition function 
(|4.17p given by 

(4.19) Z§lix,y,v,q) = Y,Q^^^ q<^^l^ 5x{x)5^{y) , 

A 

which is proven in [M] to be a tau-function of the Toda lattice hierarchy. Here k{\) is one of 
the shifted symmetric polynomials that has the explicit expression 

(4.20) k(A) = E ^i i^i - 2i + 1) , 

i 

which is essentially the Casimir eigenvalue C2 (A) [99]. It is then clear that (|4.19p can be specified 
to describe (|4.15p simply by choosing x = y = Q~^'^ (1, . . . , 1). This implies that the Yang-Mills 
partition function Zyivi (^^) ^^ ^ tau-function of the Toda lattice hierarchy. 

Once again, the g-deformed case is the one more directly connected to Donaldson-Thomas 
theory, and the formula (j4.19p applies to that case as well via a suitable specification of x = 



The genus g of the covering surface is determined in terms of the branching structure and the degree d by 
the Riemann-Hurwitz formula. 
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(xi,X2, • • • ) and y = {yi,y2, ■ ■ ■)■ The heat kernel expansion for the partition function of q- 
deformed two-dimensional Yang-Mills theory on the sphere can be written in terms of Schur 
polynomials as |100j 

(4.21) Zj:_(^),(pi) =5] B,(l,, ,...,9^-1)2 ,^^^W , 

A 

where Sa(1, q, ■ ■ ■ , q^~'^) = diuiq A is the quantum dimension of the U{N) representation associ- 
ated to A. As above, it can be interpreted as a tau-function of the Toda lattice hierarchy and it 
is also a generating function for double Hurwitz numbers. As before, due to the argument of the 
quadratic versus linear behaviour of the exponential term in ()4.2ip . it is the partition function 
of the f7(oo) chiral sector of g-deformed Yang-Mills theory that can be related to the orbifold 
partition function Z^rj.{qP, Q) (or to the instanton partition function of the K-theory version of 
the gauge theory above). 



4.5. Equivariant Gromov-Witten theory. According to the work of Gross and Taylor |10H 
I102^ I103J , two-dimensional Yang-Mills theory has a large A^ string expansion based on branched 
coverings of P^. The Hurwitz theory involved is the one that enumerates simple Hurwitz numbers 
[102^ 1103^1104) . In the formalism of double Hurwitz numbers, it corresponds to taking X2 = X3 = 
• • • = above. On the other hand, in the application of [34] to two-dimensional Yang- Mills theory 
and its g-deformation on P^ the two sets of variables x = y are taken to be equal, in order to 
describe the s| terms of the partition functions (as occurs throughout this paper). Hence only a 
particular case of (|4.19p is required, and this seems to indicate that it is more natural to consider 
two branch points at and 00, instead of just one, with the same monodromies. This is true for 
both the ordinary and the g-deformed case, changing only the branching data which depends 
on the parameter q in the latter case. 

The g-deformed gauge theory conjecturally describes a non-perturbative completion of the 
A-model topological string theory on the local toric Calabi-Yau threefold given by the total 
space of the rank two holomorphic bundle Xp = Opi (p — 2) © Opi (— p) over P-*^ (generalizing the 
resolved conifold Xi = X and the C^/Z2 orbifold resolution X2 = Y). The closed perturbative 
topological string partition function on Xp is given by |100^ I104J 



(4.22) zfj^iq,Q) = J] (-1)^1^1 QI^I q(P-^)-W/^ Wx{q) 



where 

(4,23) wM = <r''"'U^i 

is a specialization of the topological vertex [6j. In terms of the principal (7-specialization of the 
Schur functions given by (|2.7|) . the topological string partition function reads 



(4.24) zlliq,Q) = ^(-l)Pl^lQl^lg-2"W+(P-3)«W/2s;,(g*-i/2) 



A 
where we have dropped an irrelevant overall constant. 
Two distinct simplifications of ()4.24p occur for p = 1 

(4.25) Zi(g,g) = Y: (-Q)'"' Q-^^ ^?+I^U.(cz-V2)2 

A 
and for p = 3 

(4.26) Zf:^iq,Q) = Y^ (-Q)I^I g-2-(^)sA(9*-i/2)2 . 

A 
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Y 

As explained in [I04j . here Z^^{q,Q) is an equivariant partition function which describes a 
topological string theory that is generically different from the standard one. For p = 1 it gives 
the usual topological string theory on the resolved conifold X, and corresponds to a 1-Toda tau- 
function. For p = 2 the equivariant topological string theory has partition function Z^.^ {q, Q) = 
Z^Jp{q, Q)~^, consistently with the duality (j2.24p between the corresponding Donaldson-Thomas 
theories. The case p = 3 corresponds to the local P^ partition function on C'p2(— 3) — )• P^ [Ji)J]. 
In this instance there is no second Casimir term and hence no quadratic dependence on the boxes 
of the representation Aj. This seems to imply that it corresponds to a 2-Toda tau- function, in 
contrast to its non-perturbative completion, which is presumably related to the fact that the 
local P^ geometry contains a compact four-cycle. It would be interesting to better understand 
in general which topological string theories correspond to 2-Toda hierarchies, for which the 
string theory may be more naturally described using Okounkov's formalism of double Hurwitz 
numbers. 

A matrix model for the topological string partition function Z^J {q,Q) was found to leading 
orders in |104 j , and subsequently extended to all orders using the g-deformation of the Plancherel 
measure by Eynard in [105J . Explicit connections between Nekrasov partition functions and 
topological string amplitudes through the constructions of matrix models from the point of 
view of the Plancherel measure and its generalizations were found in [1061 11071 I108J . which 
rederive the results of |1U41 I1U5] from a more general perspective. These matrix models were 
recently extended to more general geometries in [28j . These matrix models for Nekrasov partition 
functions are also related to the infinite chamber limit of the Donaldson-Thomas matrix models 
found in El 



5. Vicious walkers and stochastic growth representation 

5.1. Lock-step model. In [45j, we find an interpretation of certain random unitary matrix 
model averages 

N 

ndet(/ + x.C/t)det(/ + y.[/))^^^^ 



(5,1) =n / ^n '+-<«-"'' i+^-<='"'ni<'' 



j=l •'^ 1=1 k<l 



9. _ei^f 



in terms of configurations of weighted non-intersecting lattice paths. The non-intersecting path 
model associated to (|5.1|) is the lock-step model, with a slight variation. It is explained in detail 
in [l5] and we summarize their description in what follows. On the x-axis, the allowed points 
are x = 1, 2, . . . , A^. The procedure is that each point is moved to the line y = 1, according to the 
rule that each x coordinate must either stay the same (unit weight) or increase by one (weight 
xi), always with all x coordinates remaining distinct. This procedure is repeated a total of N 
times, with each right diagonal segment at step j weighted by Xj as x coordinates are moved to 
the line y = j- Notice the difference with the original lock-step model, where at each time all 
particles move either to their right or to their left with equal probability. After step N, perform 
another N steps, but now with the segments either vertical (unit weight) or left diagonal (weight 
y2N+i-j at step N + j). As usual, we have the conditions that the segments do not intersect 
and, in addition, are further constrained to return at the line y = 2N to the same initial x 
coordinates. 
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We thus have, as was the case in the Brownian motion description of Chern-Simons partition 
functions [lOj, a reunion condition on the walkers. The result of [l5] is that the generating func- 
tion for this process is given by (|5.ip . As we have seen in Section 13.21 in the limit A^ — )• oo with 
the weights of the process taken to be Xi = yi = g*~^'^, this is the Donaldson-Thomas partition 
function Zj^r^ (q). Since N = oo there are infinitely many walkers. Likewise, the reduced orbifold 
Donaldson-Thomas partition function Z^rj,{q, Q) is the generating function for this process with 
weights Xi = yi = \/Q ("9)*"^ • More general geometries correspond to higher-dimensional ver- 
sions of this process, with the walkers constrained to move independently on coordinate planes. 
In this picture, wall-crossing phenomena appear as the creation or destruction of coordinate 
planes, and hence particles, and changes in the weightings of walkers. 

The celebrated Robinson-Schensted-Knuth correspondence |51j establishes a bijection between 
weighted N x N non-negative integer matrices with entries aij weighted by {xi yj)"'"'^ and pairs of 
weighted semi-standard tableaux of content N . This implies that (j2.2p is the generating function 
for the weighted matrices, and hence the MacMahon function M(q) is the generating function 
for infinite-dimensional matrices with entries aij weighted by (^g(*~^/2Hi-i/2)a-ij"j _ Moreover, in 
it is found that the function 

N 

(5.2) J] (l-x^xj)-^ 

is the generating function for 2A^ x 2A^ matrices [aij] which are invariant under reflections about 
the entry {N + ^,N + ^), i.e. aij = a2N+i~i 2N+i-j- If we specify Xi = q^~^/^ and take the limit 
N —7- oo, the generating function is the square of the MacMahon function, a factor that appears 
in the Donaldson-Thomas partition functions of the conifold and the C^/Z2 orbifold. 

Interestingly enough, instead of considering a symmetry constraint on the matrices, one can 
impose constraints on the entries. In particular, if we constrain them to be or 1, with the 
entries aij weighted by (xiyj)"''^ , then the generating function is the right-hand side of the 
dual Cauchy identity (|2.19p [45]. Recall that this leads to the Schur expansion (|2.20p of the 
Donaldson-Thomas partition function on the resolved conifold. Thus the reduced partition 
function Zf^rj, (g, Q) can be interpreted as the generating function of such matrices, with weights 
Xi = —\/Q {.—<if^ and yj = —VQi—q)-'^^''^- These alternative enumerative descriptions in 
terms of infinite non-negative integer matrices are intriguing, in light of the fact that in general 
the Donaldson-Thomas partition function Z^rj, of a threefold X is a generating function which 
counts ideal sheaves on X. It would be interesting to formulate a more direct connection between 
these infinite-dimensional random matrix theories and the geometric counting problems. 

5.2. Corner growth model. It is also possible to relate the Donaldson-Thomas partition 
functions to the distributional limit of the corner growth (or last passage) model with geometric 
weights [18]. Let u {i,j), {i,j) S ^^ be independent geometric random variables and define 

(5.3) G{M,N)= max V u;{i,j) , 

where the maximum is taken over all up/right paths vr from (1, 1) to (M, N); this is called the 
corner growth model [18]. It is shown in |19) that if w (i,j) are independent geometric random 
variables with probability distribution of the form 

(5.4) V [oj (j, k)=m]={l- xj yk) {x, y^T , 
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then the probabihty that G{M, N) is smaller than a certain value can be written in terms of the 
Schur measure as 

n 

(5.5) V[G{M,N)<t] = (^l{{l-xjyk)) ^ Sa (x) Sa (y) . 

j,k=l \:\i<t 

Since it is given in terms of the Schur measure, the choice of parameters Xj = q^"^'"^ and y^ = 
(jrfc-i/2 niakes the normalization constant of the process in the limit n — )• oo (the normalization 
constant for t — )• oo) equal to the Donaldson-Thomas partition function for C^ in ()2.ip . 

As in the case of the lock-step model, since the expression (j5.5p is very general as it in- 
volves generic coefficients Xj and y^, it also leads to the orbifold partition function (j'2.'25p by 
again choosing Xj = \/Q {^—qf~ and y^ = ^/Q {—q) ~ . To describe the conifold partition 
function (j2.20p . we need to construct the dual process and hence modify ()5.4p to 

(5.6) P^ [to (i, k)=m]={l + Xj VkY^ Sfe(x,- y^T , 

where Hf, is the endomorphism of the ring of symmetric polynomials given by H^ (ca) = f)A for 
every Young diagram A, with e\ the elementary symmetric polynomials and f)A the homogeneous 
symmetric polynomials [51]. Then 

n 

(5.7) P^ [G(Af , iV) < t] = ( n (1 + ^J y^)'' ) E ^A (^) ^A' (y) > 

j,k=l \:\i<t 

and the choice Xj = —\/Q {—qy~ and y^ = —\/Q {—q) " makes the normalization constant 
of the process equal to the partition function (j2.20p in the limit n — )• oo. More general geometries 
correspond to multiple copies of this process involving independent random variables. Wall- 
crossing in this picture is the creation or destruction of independent random variables and 
changes of weightings. 

A particular case of this stochastic growth model can be understood as a generalization of 
the stochastic process underlying the longest increasing subsequence problem. The distribution 
of the Poissonized version of the random variable L(a) describing the longest increasing subse- 
quence in a random permutation is given by the Gross- Witten model |109j . The choice of the 
distribution (|5.4p is a generalization of the original model introduced in [lU] which has Xj = x, 
yk = i for all j,k. If one takes x = a/N"^, then G{N,N) converges in distribution to L{a) as 
A^ — )• oo, and so one can view G{N, N) as a generalization of the random variable L(a). 

This model is intimately related to other growth models, as discussed in detail in |18[I19). A 
growth model is a stochastic evolution for a height function h{x,t), with x denoting space and 
t denoting time. An admissible height function has to satisfy h{x + l,t) — h{x,t) = it 1 for all 
t. In particular, the discrete polynuclear growth model is a local random growth model defined 
inductively by 

(5.8) h{x, t + 1) = max {h{x - 1, t) , h{x, t) , h{x + 1, t)) + uj{x, t + 1) 

with {x,t) G Z X N and h{x,0) = 0, where u}{x,t) are independent random variables. One can 
think of h{x, t) as the height above x at time t, so that the map x i— t- h{x, t) describes an interface 
evolving in time. The special case where uj{x, t) = if t — x is even or if |x| > t, and 

(5.9) w{i,j) =uj{i-j,i+j -1) 

for {i,j) G Z^ are independent geometric random variables with probability distribution ()5.4p . 
yields G{i,j) = h{i—j, i+j—1). This growth model is expected to fall in the Kardar-Parisi-Zhang 
universality class. 
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It is also possible to demonstrate the equivalence with other systems, such as random tilings of 
Aztec diamonds (which is related to a dimer model), and non-intersecting walks on a graph |18t 
[19]. The latter correspondence, in the case of the distribution (j5.4p that leads to Donaldson- 
Thomas theory, is known in detail [TS] but its description is rather lengthy. It seems that the 
lock-step model |45] described above should be directly related to this version of the corner 
growth model since, as we have seen above, the probability ()5.5p without the normalization is 
the generating function of that vicious walkers model. 
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